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- z2— concentration of emitted nitrogen oxides Cyox [ppm]; 
- z3— concentration of emitted hydrocarbons Cyc [ppm]; 
- Z4—mean indicated pressure pm; [MPa]; 
- Zs — indicated power N; [kW] 
3. Set of values C, which can appear during the test, were accepted as constant (fuel 
type). 
4. Set of interfering values H, which can be modified (ex. ambient temperature, 
atmospheric pressure, relative ambient air humidity). 


Due to the low number of input values in the research, there was applied a complete 
three-valued 3° plan with 3 independent variables and 27 measurements carried out in a single 
measurement block. It is presented in Table 2. The scope of three fuel temperature values 
covered two extreme values, imposed by physical and chemical parameters of fuel as well as 
by the safety conditions during increasing the fuel temperature. The third value was the mean 
value of the two extreme temperature values. Other input values, such as engine speed and 
torque, were limited by the engine construction parameters. The extreme values and their 
arithmetical mean were accepted as in case of fuel temperature. 


Table 2. Plan of the experiment 3** - fractal complete three-valued plan I block, 27 systems and received 
results of measurements 


Input data EANA data 
measureme - 
nt system c = co = = HC Pmi Ni 
El p: ni [ppm] | [ppm] | [ppm] | [MPa] | [kW] 


937.8 191.7 137.4 0.410 9.00 
1079.3 208.9 173.9 0.577 12.67 
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2.2. Statistical and content-related analysis of the test results 


442.9 


22 
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25 
26 
37 


126.5 
53.7 
80.9 
93.5 
112.6 
159.2 
149.1 


0.801 17.57 


Registered parameters were statistically analyzed at confidence interval a = 0.05. After 
performing analysis of match measures of the significant statistic coefficients, there was 
selected an interaction L-L (linearly-linear) that had the most beneficial values of the quotient 
of regression coefficient and mean error of estimation the regression coefficient t = bi/Spi. 


For the accepted model there were plotted some graphs presenting the estimation of 
match effects and the influence of the respective input values on a considerate output value. 
The graph of standardised effects is an effective tool presenting which factors have the 
greatest influence on the output value. As an example, Fig. 4 shows the influence of fuel 


temperature on the concentration of HC in exhaust 


gases. 
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Fig 4. The graph of the standardised effects of HC concentration in an outlet manifold of the engine tested 
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Wartość aproksymowana 
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Fig 5. Dependence of the measured and approximated values of HC concentration in an outlet manifold of the 
engine tested — total of residues MS = 423,2653 


As it results from Fig 4, the greatest influence on the HC concentration have: fuel 
temperature, followed by torque and the engine speed. 

Torque in the first place and engine speed in the second one, are the most significant 
parameters that affect the NOx concentration. Fuel temperature does not play an important 
role in the NO, formation in the combustion process. The greatest influence on mean 
indicated pressure has torque and, in the definitely lower extends — the fuel temperature. The 
greatest influence on mean indicated pressure and indicated power has in order: torque, 
engine crankshaft speed and fuel temperature. Fig 5. illustrates the dependence of the 
measured and approximated values of HC concentration in a outlet manifold of the engine 
tested. 

Engine torque has the most significant influence on CO concentration in a outlet 
manifold. At torque boost, the dose of fuel increases, air excess coefficient 2 decreases and it 
results in the complete and incomplete combustion. Increase in torque 7 at constant air excess 
coefficient A goes with the decrease of CO concentration. The decrease in torque 7 causes 
decrease in combustion temperature, which also involves the decrease in CO concentration. 

Decrease in £,, results in deterioration of the quality of fuel pulverization and injected 
fuel droplets diameters worsen the condition of fuel evaporation, which is caused also by the 
temperature decrease in the centre of fuel injection. Then, self-ignition delay increases, which 
causes interferences in the combustion course and as a result increases the amount of products 
of incomplete combustion, including the increase of CO in the engine exhaust gases. 

The amount of toxic compounds in exhaust gases in the test engine increases along 
with engine load as a result of decrease in A and increase in combustion temperature. When 
torque T (from T = 25 N-m to T= 125 N-m), at n = 800 [rev/min] increases, CO concentration 
increases as well — Fig 6. 
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The greatest influence on the increase in NO, concentration in four cycle engine has 
the torque (7). It results from the fact that two factors: temperature and amount of oxygen, 
playing a significant role in NO, formation, affect one another. Along with the increase in 
torque T, the value of A decreases, which results in temperature increase at the simultaneous 
decrease in oxygen concentration Along with the increase of toque 7 (from T = 25 to T= 125 
[N-m]) at n = 1000 rev/min, the air excess coefficient A decreases., which results in the 
increase in oxygen concentration in exhaust gases. 

In general, an increase in temperature of injected fuel leads to decrease in NO; 
concentration in exhaust gases, which is the most noticeable in the temperature range 
50+100°C. Thermal activation of fuel is the most efficient at low loads, because along with 
the increase in engine load, the effectiveness of this procedure radically changes and leads to 
a completely inverse situation, which is illustrated in Fig 7 
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As it was mentioned above, the principal influence on the increase in HC 
concentration in a outlet manifold has the temperature, what was verified in the statistical 
analysis. It is shown in Fig. 8. 
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It was possible to receive the decrease in uncombusted hydrocarbons along with the 
increase in fpa; (even to approx. 55% of the values received at tpa = 20°C). However, it 
depends on engine work conditions: an increase in both its load and torque powers that 
beneficial effect. 

The increase in HC concentration in exhaust gases, as a result of increase in torque T , 
at a given engine speed n, results from the decease in air excess coefficient A. It is assumed, 
that in some areas of the combustion chamber there can appear local decrease in å to value 
(A<<1), which has the significant influence on increase in the number of uncombusted fuel 
particles. As a result, it leads to increase in HC concentration in exhaust gases. It is noticeable 
that HC, similarly to CO, mainly forms as a result of heterogeneous air-fuel mixture. This 
heterogeneity, in high extend depends on the structure of combustion chamber and its system 
of air filling and fuel feeding. 
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RAS: Fig 9. Change in mean indicated pressure pmi as a 
n = 1200 [obr/min] E 8i 
a i function of tha = (20+150)°C at a given engine speed 
praet iasa and variable torque T 
0,7 
= a 
& 
2 0,6 


Torque has the main influence on the increase in mean indicated pressure in engine 
cylinder. Second in order is the fuel temperature ż,„ (Fig 9). 

Equation of the regression of CO, NO, and HC concentration, mean indicated pressure 
Pmi and indicated power N; in the range of tu = (20 + 150)°C, n = (800 + 1200) rev/min and T 
= (25 + 125) N-m are the following: 


[Cco] = + 252.8545 — 19.8306, + 0.0703 tyai + 0.8471n — 0.0003» + 15.01777 
+ 0.0658T”+ 0.0086ż,am — 0.073 6taiT — 0.0068nT 


[Cox] = — 228.7217 + 0.911 6f,ai— 0.0033 tpai + 0.6561n — 0.00037" + 3.54227 
+ 0.00457" — 0.0006ż,am + 0.0053 t,aiT — 0.0027nT 


[Cuc] = 479.2674 — 2.7329tpar + 0.0056tpai — 0.7529n + 0.00040" + 2.1907T — 0.00667" 
+ 0.0014 tain — 0.0011 tpaiT — 0.0005nT 
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[Pmi] = 0.3246 — 0.0014tpa + 0.0000tpa” — 0.00017 + 0.00007? + 0.00507 — 0.00007” 
+ 0.00004, + 0.00004,,77 — 0.0000nT 


[Ni] = — 6.3612 — 0.0278tpai — 0.0000%, + 0.01717 — 0.0000n” + 0.08447 — 0.00017" 
+ 0.0001 tpa — 0.0001 tyaiT’ + 0.00007 


Dependence of NO, concentration on the values of input values (including fuel temperature) 


is presented in Fig 10, whereas the dependence of mean indicated pressure on these values is 
presented in Fig 11. 
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Fig 11. Mean indicated pressure pm; as a function of tpa and T 


The highest increment of mean indicated pressure value (28%) for heated fuel is at n 
=1200 rev/min and T = 25 N-m. Low increments of this pressure are registered at low engine 


speed and high load. The increment is low and it reaches only (6%) in regard to the same load 
at unheated fuel. 
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3. Conclusions 


The main aim of the use of thermal activation of fuel was to estimate energy and 
ecological effects of heating fuel in the feed system of vessels combustion engine. 
Considering the influence of temperature of injected fuel on analysed engine performance 
index, it is rational to search its optimal value. 

Mentioned above factors must be considered up to the arbitrary decisions while 
making an attempt of finding the optimal values £,,,. Also engine work conditions must be 
taken under consideration while comparing the results for particular values tpa 

On the basis of test results it can be stated that: 

1) Along with the increase in engine load, there is a decrease in fuel temperatures, 
required for the most beneficial effects of heating the fuel. 

2) Optimal fuel temperature in the range of (80+120)°C apparently decreases the power 
input necessary to its receiving in proportion to the input required at £,„ = 150°C (they 
comprise at average (80+90)% of the input at heating up to 150 °C). The range of 
temperature given above constitutes also the most beneficial range in proportion to the 
decrease in toxic compounds emission. 

Feeding the engine by heated fuel brings on generally beneficial results. During 
proceeded tests there were obtained the following decreases in: 

- HC concentration for about 44% (100 ppm) 

- Carbon monoxide CO concentration for about 40% (1200 ppm) 

- NO,, concentration for about 22% (50 ppm), but only at low load, because at high load 
NO,, concentration increases for about 20% (50 ppm). 

Considering energy parameters of the engine fed by heated Diesel fuel in comparison 
with cool fuel, there was remarked the increment of: 

- mean indicated pressure for about 14% (0.1 MPa) 

- indicated power N; for about 13% (3 kW). 

Comparing research conducted on the test stand in LEUO with tests carried out by I. 
Pielecha from Poznan University of Technology it is possible to state that the obtained results 
do not differ much, although in the research of Poznan University of Technology the fuel was 
heated up to 220°C, while in the research of Naval Academy of Gdynia it was heated up to 
150°C. 
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Abstract 


The following article presents the possibility of enter variables aggregation in neuron models required 
power for the seagoing vessel by means of dimensional analysis. Such aggregation simplifies significantly the 
model, on the basis of which, the required power for ship propulsion is defined. The algebraical diagram of 
dimensional analysis constructed by S. Drobot is a good instrument for this goal. What is more, such diagram 
allows to control, in respect of mathematics, the correctness of conclusion rules used in neuron models. 
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1. Introduction 


Power for ship propulsion is selected on the basis of hull resistance analysis in different 
conditions of sailing. At present it is possible to observe more and more the use of artificial 
neuron networks in selecting power for ship propulsion [1,2]. 

In diagram 1 we can see an exemplary structure of neuron network allowing to determine 


power for ship propulsion in different conditions of sailing. 
enter 


Fig.1 Structure of neuron network determining required power for ship propulsion in different conditions of 
sailing: M- torque of the propeller tail-end, n- revolution of the propeller, t- time of change of respective 
parameters, v- speed of ship, I- inertia moment of whirling propulsion elements, D-mass of ship 
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Therefore ship models are created using activity of neuron networks. Such networks have 
a complex structure putting into practice power selection for ship propulsion on this basis in a 
very wide range of variable parameters. The power of ship propulsion depends on dimensional 
quantities given to enter part of network which is put into activities respective practice. These 
activities involve checking premises of rules enclosed in respective bases of measurement 
data. 

Conclusions of rules concerning power selection for ship propulsion are defined by 
neuron network. The number of power selection rules for ship propulsion increases violently 
according to the increase of the number of entrances [3]. 

In order to reduce the number of enter signals being dimensional quantities it is possible to 
use their aggregation by means of dimensional analysis. Aggregation involves the replacement 
of enter signals by the signal which is a respectively selected product, and their functional 
combination. 


enter 


exit N 


bias 
m 


Fig.2 Structure of neuron network determining ship propulsion power followed by their aggregation of enter 
quantities: p-nondimensional number of time,, p,-nondimensional number of speed, py—nondimensional 


number o propeller torque, y = -nondimensional number of power requirement 


nel 


To achieve good results from aggregation it is primary important to select proper 
dimensional function combining enter variables. At the same time full use is being made in 
physical relations between power required for ship propulsion and ship resistance. Ship’s hull 
resistances in different conditions of sailing are enter quantities into neuron network. In fig.2 
the diagram of neuron structure network determining propulsion power of the ship after 
aggregation of its enter quantities, has been presented. 
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2. Aggregation of dimensional quantities describing requirement of power for ship 
propulsion 


Power requirement for ship propulsion can be described by means of the following 
dimensional function [4,5] : 


N=GM, n, t,v,l, D) (1) 
where: 


N — required power for ship propulsion in [ kg m” s” ], 

M — torque at the tail-end of the propeller in [ kg m” s” ], 

n- rotational speed of the propeller in [ s" ], 

t — time of change of respective parameters in [ s ], 

v — speed of the ship in [m s" ], 

I— inertia moment of whirling mass of the ship propulsion system in [ kg m” ], 
D — mass of the ship in | kg |]. 


The function defined by the formula (1) is a dimensional function in a dimensional space, 
whose arguments are its quantities. Such function cannot be treated as ordinary numerical 
function in the set of real numbers. Dimensional function (1) fulfils the conditions of 
invariance and dimensional homogeneity (6). Such conditions do not limit the form of 
numerical function where arguments as well as the value of the function are nondimensional 
quantities. 

Six argument of dimensional function of required power for ship propulsion (1) have 
respective dimensions in an adopted system of measurement units SI. The matrix of exponents 
of the dimensions in this case, has the form: 


Kg m S 
M 1 2 -2 
N 0 0 -1 
T 0 0 1 
V 0 1 -1 
I 1 2 0 
D 1 0 0 


The above matrix is of the third order, which means that among arguments of dimensional 
function three arguments are dimensionally independent. 

Dimensional analysis does not provide any method which could make it possible to 
ascertain which quantities can be chosen as dimensionally independent. Generally, if six 
quantities define propulsive power of the ship and three quantities are dimensionally 
independent so the number of possibilities of their selection equals three elements 
combination out of the six elements set resulting in the number of twenty. Those combinations 
are, by way of example, the following quantities dimensionally independent of the remaining 
ones:( n, I, D), (t, v, D ), (n, I, t ), ... etc. After selection of arguments dimensionally 
independent of the remaining ones, we can, on the grounds of Buckingham 's theorem ( tw II ) 
write down the dimensional function (1) in the following forms: 
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N= f (Ps pu) në: P: D° (2) 
where: 


f — numerical function defined in the domain of real numbers, 

P+, Py, PM — nondimensional element of respective dimensional quantities, 
a,b,c,d — real numbers, 

- the remaining denotations as in the formula (1). 


Exemplary relations between dimensional quantities dependent: t, v , M on independent 
quantities: n, I, D have the following form: 


t=: n P-D v= ny If e DY M=ogy'nż' 1) - D’ (3) 
where: 


a, 01,02,8, Bi, B2 Y, Yi, Y2 — real numbers, 
- the remaining denotations like in the formulas (1) and (2). 


Real numbers being the exponents of dimensional quantities are determined by 
comparison of dimensions of quantities taking part in relations (3) and (2) in the following 
way: 


[s]=[s" ]*: [kgm]; [kg T > Pep WO 
[ms ]= [s [i [ken [ksh + wsi; b=: n=- 
[kgmðs?] = [s7]: [kem] [ke] > œs} a= n= 0 
[kgms] = [st [kgn [kg >  a=3; b=; c= 0. 


Calculated values of numerical exponents taking part in equations (2) and (3) allow us to 
write down the equation (2) in the following way: 


T A (4) 


where: 

- denotations like in formulas (1). 

The formula (1) proceeding in the same way as above we can obtain different numerical forms 
of dimensional functions (1) while selecting other quantities as dimensionally independent. 


So, selecting the quantities: t, v, D as dimensionally independent of the remaining ones, we 
can obtainces: 
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ak =f Le BE. (5) 


where: 
- denotations like in formulas (1) and (2). 


Forms of numerical dimensional function (4) and (5) with the same exactness define 
power required for ship propulsion. There are as many such forms as there are possibilities of 
selecting quantities dimensionally independent, correct in respect of mathematical precision. 

In the case under investigation there are twenty of them. That means that during selecting 
numerical form of dimensional function we should concentrate on the way which makes it 
possible to carry out an experiment and train of the network on the set of measuring data. 
Definition exactness of such functions depends on evaluation of constant coefficients during 
experimental research by means of regression. 

Form of above named numerical function define base for perennial day and night data 
recording the state parameters of vessel movement in task: ,, Realization of systematic 
investigations of ship energy propulsion demand ” government department problem [8]. An 
exponential form of numerical function in which constant coefficients determine measuring 
data by means of least square method is selected. 

Demand of power for ship propulsion determined by means of neuron model taking 
advantage of dimensional function (1) is in possession of six enter quantities and the same 
number of conclusion rules. Neuron model of ship propulsion system contains bases of rules 
making use of numerical form of dimensional functions which can be easily interpreted and 
modified. 

Application of algebraical scheme of dimensional analysis constructed by S. Drobot 
allows to aggregate conclusion rule defined by the function (1) to form (4) or (5) in which 
three enter signals take place. So we can see that by means of dimensional analysis it is 
possible to reduce in an essential way the number of entrances into the neuron model of ship 
propulsion system and thus considerably simplify its structure. In the model under 
investigation there are only three enter variables instead of six. 

At the time of identification and also estimation of enter signal in the mode on-line of 
such models there is the Reed preliminary calculation of the value of aggregated variables. 


3. Conclusions 


Determining power required by the ship means of neuron models allows to train the 
network in the set of measuring data. Such models are represented by the base of rules for 
which algebraical scheme of analysis can be used. 

It allows to aggregate the model in different ways depending on the need aggregation of 
neuron model by means of dimensional analysis leads to a considerable reduction of the 
number of its entries. 

In this way the structure of neuron networks is simplified and used in aggregated neuron 
models. 

Precision of numerical function estimation (5) used in aggregated neuron models equals 
the exactness of determination of its parameters by way of experiment. 

Besides the dimensional analysis allows to ascertain the structure of conclusion rule, and 
being numerical function, it is properly constructed in respect of mathematical correctness. 
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Abstract 


This paper presents a concept of method intended for the assessing of operation of engine which works under 
partial load. A valuation ( quantitative) approach to the operation interpreted as a physical quantity was applied to 
determine existing margins of operational parameters of engine as well as additional costs associated with 
degradation of its technical state, born by engine’s user. 
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Introduction 


From the point of view of operation aims generated by operational system in the form of 
transport tasks for ship considered as a technical system, ship propulsion system is intended for 
realization of assumed sailing speeds of the ship within a given period and determined range of 
changeability of external conditions. 

Hence the functioning of the ship in steady motion with speed v is determined by description 
of a system of forces acting on ship’s hull, that can be analytically represented ( on assumption of 
neglecting the gravity and buoyancy forces) as follows [3, 10]: 


T=R,+AT (1) 


gdzie: 

T — thrust of the propeller, 
AT — thrust deduction, 

Rr — total hull resistance. 


The interdependence of ship resistance and effective power of main propulsion engine for 
steady sailing speed is described by the following equation [3, 10]: 
yo 2) 
Niw ` Ns 
where: 
v — speed of ship, 
Nw — efficiency of shaftline, 
1s — propeller efficiency. 
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Knowledge of the quantities appearing in Eq.(2) makes it possible to properly select all 
elements of propulsion system in the designing stage [9, 10]. 
Considering the sailing speed to be an independent variable which definitely determines two main 
factors associated with realization of objective function, namely : 

e duration time of realization of transport task, which definitely influences costs of its 

realization and potential profits to ship’s owner, 

e possibility of maintaining an assumed ship course in heavy weather conditions, which first 

of all influences occurrence probability of an emergency situation, 

one should observe that during ship operation all remaining quantities undergo changes due to 
degradation processes which involve evolution of the process of technical state changes towards 
technical unserviceability states. It obviously influences serviceability of system’s elements and 
,consequently, of the entire system. 

On assumption that main propulsion system of ship has been correctly selected, the hull 
resistance characteristics Rr = f(v) constitute the factor which generates demand for effective 
power developed by propulsion engine for reaching a given sailing speed of the ship [4]. The 
characteristics determine, for assumed conditions, value of developed effective power of engine 
neccessary to reach a demanded sailing speed. If a relatively short period is taken into account then 
it can be assumed , without any large error, that the representation is invariable. 

In practice, an important change of the relation Rr = f(v), resulting from influence of the 
environment and physical ageing processes ( especially of wear, e.g. corrosion ), occurs , whose 
external symptom is a rise of value of the force Rr at given values of the speed v in comparable 
sailing conditions. Long time intervals between successive ship hull docking surveys ( possible 
repairs) are of special importance because such intervals , depending on a type ( class) of ship and 
given classification rules, vary from 3 to 5 years [e.g. 12]. 

In selecting the main engine in the stage of ship designing the above desrcibed processes are 
taken obviously into account, that is manifested by skillfully selected design characteristics of ship 
hull resistance as well as by application of the so-called sea margin during determination of the 
main engine contractual power Nx [9, 10], which makes that the main engine usually has certain 
power surplus as compared with assumed design values of hull resistance. It means that in the case 
of typical transport task for a given ship its main engine will be used in a load state lower than 
rated one [1]. 

However, the existing power margin which guarantees ship safety, is systematically decreasing 
due to the above described phenomena and the worsening of operational characteristics of engine 
itself. 

Hence, it seems rational to elaborate a method for determining the existing margin of 
operational engine parameters associated with degradation of its technical state, and the additional 
costs born by its user. 


Influence of change of technical state of engine on its operation 


Worth mentioning, in the light of the above described power surplus of main propulsion 
engine, that for the engine running under partial loads the process of decreasing the available 
output power will develop in two phases: 

e in the first phase hourly fuel oil consumption will rise solely ( at a relatively constant value 

of developed torque ), and in consequence operational costs will be also higher; 

e in the second phase a limitation of value of effective power developed by the engine will 

occur because of constructional constraints and lack of possibility to increase fuel charge. 

The described phenomena result from control action of fuel apparatus which will increase, 
within a given range of values, the instantaneous fuel charge gp (gy ” - fuel charge for i% load of 
engine being in the technical serviceability state under assumption that the maximum engine load 
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amounts to 110% of its rated value, i <//0) up to the instant of reaching its maximum value 
Gpmax. Every successive decrease of value of engine total efficiency will result in a recordable 
decrease of its torque Mo. 

If to assume the engine partial load to be constant ( e.g. for 9RT-flex60C-B engine : 85% 
engine load, the contractual output power Nx = 80% Nri = 17370 kW, the contractual engine speed 
nx = 90% ng; = 102,6 rpm [13]) the phenomenon in question can be presented in the form of the 
diagram (Fig. 1) of the following interpretation : 

e the time-dependent drop of the engine's total efficiency (in the case in question by about 
9%) results , in the first phase, mainly in an increased hourly fuel oil consumption ( 
specific fuel oil consumption). It can be represented as the occurrence of the successive 
recordable events F which consist in increasing the fuel charge go” by the value Ag, at a 
relatively constant , relevant to a given engine load state, value of the torque Mo. This way, 
increased operational costs are generated rather without imposing limitations on the ship 
motion parameters described by Eqs. (1) and (2). 

e the degradation processes progressing along with time of further use of engine result in 
occurrence of the recordable events U which consist in decreasing values of the engine 
torque M, at constant fuel oil consumption (at gy = Gy max). Further long-lasting use of the 
engine results in the significant worsening of its characteristics which definitely impose 
limitations on ship motion at an assumed speed or course. In heavy weather conditions 
such situation will obviously form an important circumstance for occurrence of state of 
danger to ship safety. 
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Fig. I Changes of engine fuel charge and torque , which result in changes of total efficiency of the engine 
running under partial load 


Assessment of available power margin and additional operational costs of engine 


The number of repetitions of the event F, N4gp , within the time interval (0, t) is a random 
variable of non-negative integral values. The dependence of the random variable on time 


N(t):t > 0} 


constitutes the stochastic process t . On the assumption on its stationarity, lack of 
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consequences and flow singularity [5] the Poisson homogeneous process [2] can be applied to 
description of the process of increasing the fuel charge gp” (6) resulting from the decreasing of the 
total engine efficiency Ne (in maintained engine load conditions ). 

Hence, in the case if the number of the events F „Ng, , occur up to the instant ¢, the total 
increase of the fuel charge gy! % up to the instant, can be expressed as follows: 


AG, = Ag Nagp (3) 
where: 
AG,- total increase of the fuel charge gp after occurrence of N4g, number of the events F, 


Ag, — elementary recordable fuel increment by which the fuel charge gp” increases , whereas 
the random variable N4gp has the following distribution [1]: 


P(N sp =k)= a exp(-A,1); k=1,2,...,n (4) 


where: 

Aę — a constant interpreted as the intensity of occurrence of the event F ( i.e. increase of the fuel 

charge gp” by Agp value). 

Assuming that the unit fuel charge increase by Ag, value generates the additional unit cost p 
(realization of the same task at an increased fuel oil consumption) one can determine its total value 
®, up to the instant ¢ , as follows: 

D = P : Agp’ Nagy (5) 

For the engine load smaller than maximum, the summary fuel charge increase AG, ma and 


pmax 
occurrence of certain number , N'agp , of the events F affecting solely the increasing of the unit 
fuel oil consumption ge and hourly fuel oil consumption By , can be determined by the following 
relation : 


A G + gp á =; Agp: N’ 4gp a or ina GAż (6) 


pmax 


where: 
g” — fuel charge relevant to a given engine load i % (under assumption that the maximum 
engine load is equal to 110% of its rated load, i.e. i</10) 
G max ~ maximum fuel charge possible to be realized by injection apparatus within existing 


pmax 
constructional and control limitations. 

Occurrence of the successive event F (Nagp = N’ agp + 1) for the current engine load i % will be 
associated with occurrence of certain limitation in developing the demanded power of the engine , 
and in consequence limitations for motion of the entire ship, will occur. 

As the effective output power of engine can be expressed as follows [11]: 

Te Wa Sp, 'N (7) 
T 


N, = 


where: 
Ne — effective efficiency of engine, 
wa — fuel lower calorific value, 
gr - temporary fuel charge (for i% engine load), 
n — rotational speed of engine, 
1 - a coefficient which takes into account number of strokes realized during one engine 
working cycle, 
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as well as on assumption that the changing process of engine technical state is continuous over 
time and states, the situation can be illustrated in the way shown in Fig. 2. 
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Fig. 2 Change of engine operational characteristics, which results from degradation of its technical state. ty) — 
demanded time of realization of the task 1, D,, — demanded value of operation for realization of the task 1, ty2— 
demanded time of realization of the task 2, tm; — possible time of realization of the task 1 in the situation A (described 
in the text), tmz — possible time of realization of the task 2 in the situation B (described in the text), tn3 — possible time 
of realization of the task 2 in the situation A (described in the text), D2 — demanded value of operation for 
realization of the task 2, Dm — possible operation of engine in the situation 2 , D’ — possible operation of engine in 


the situation A. 


By introducing the notions of the demanded operation Dw and the possible operation Dm , 
whose detail interpretation can be found in [6, 7, 8], the phenomena graphically presented in Fig. 2 
can be highlighted as follows: 


in the case of a hypothetical lack of possibility of increasing the fuel charge ( the situation 
A ) a drop of total engine efficiency will result in a sudden drop of its effective power — the 
line N”(t) — that , apart from producing a situation hazardous to ship safety , during further 
degradation of technical state of engine , will cause extending the realization time of 
transport task from the instant ty to tm, 

if at engine partial load it is possible to increase the fuel charge within the interval of 


values (gr, G ) (the situation B) , and up to the instant of reaching its maximum value 


pmax 
(the point A, instant t’) no noticeable drop of power occurs- the line N’(t), then only an 
additional cost resulting from increased fuel oil consumption will be generated. Hence it 
can be stated that if a transport task lasts shorter than up to the instant t’ any lengthening of 
its realization time will not occur. 

begining from the instant t’ ( corresponding with the point A) a limitation of the possible 
engine operation Dm [6, 7, 8] occurs , that is associated on one hand with the lengthening of 
realization period of transport task (or lack of possibility of its realization during a given 
period) on the other hand — with increased costs of engine's use. 
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By using the relation (4) and (6) the number of the events F , k’, corresponding with that of 
N’ agp Can be expressed as follows : 


_ oi% 
k'= Cm Sp (8) 
Ag 3 
and its occurrence probability as : 
RTA 
A „| Ag, | 
PIN = = BEAL) 21) 0) 


i% 
CA a Ep 
Ag, 


Practically possible determination of stochastic process parameters 


From practical point of view the crucial problem is to determine Ag, and Ay values. This will be 

possible if two complementary conditions are satisfied: 

e to have access to results of operational investigations carried out with application of 
standard control measuring instruments and a system for diagnosing the process in 
question by recording and analyzing changes in values of fuel flow rate within time 
intervals when engine load is relatively constant; 

e to analyze technical documentation of the engine and carry out simulation investigations in 
order to elaborate a mathematical model of influence of the occurring events F on 
operational fuel system parameters. 

In the presented case Ag, can be much easier determined , because of its small value, under the 

assumption that : 


ga PT (10) 
n 

where: 

By — unit fuel oil consumption of engine, 

n — rotational speed of engine, 

1 - a coefficient which takes into account number of strokes during one engine working 

cycle. 

The relation (10) at assumed class of accuracy (e.g. 0,5) of commonly applied flow meters 

makes it possible to determine Ag, value with an accuracy sufficient for practical purposes. 

Much greater difficulties are associated with determination of A; value, which results first of all 
from lack of empirical research in this area of engine functioning. However such value can be 
determined, in a way sufficiently accurate from practical point of view, by means of the following 
procedure: 

e on the assumption that a typical ( most frequently occurring ) load state of engine is x% [1], 

the difference between its maximum load and the above mentioned state determines the 


changeability range of instantaneous fuel charge (Bes Cm)? hence on the basis of Eqs. 
(3) and (6) it yields : 
AG, =G y (11) 


pmax pmax Sp 


e the period TaGpmax, during which AG „max value is reached, should be assessed on the basis 


pmax 
of technical documentation of engine and its producer's recommendations for overhauls 
and repairs, as results of relevant empirical research are usually lacking; 

e taking all the above into account one can state that within the considered period about 


162 


W. ie r 12 

Ag, Ag, ( ) 

number of the events F will occur, hence : 

x% 
G pmax = Sp 

N' A 

ję a (13) 
TG Pmax TG Pmax 


Finally , making use of the above presented procedure and basing on the relation (8) one can 
determine at first the number of the events F, k’ , and the probability of their occurrence: 


ze] 
Grieg z 5 * Asp 
Ag p P 
TAG Gia Aj 
^g, 14 
P(N sep k') > exp t > (14) 
Gz z Ep pom 
Ag > 


and, the application of the relation (14) is this much practical that in the case of lacking results 
of empirical research one can approximately (but with a sufficient accuracy) determine 
P(N hep = k')value , basing solely on engine’s technical documentation. 

For example, for the above mentioned 9RT — flex60C — B engine (at the contractual output 
power N, = 80% Nr; = 17370 kW, the contractual engine speed n, = 90% npg; = 102,6 rpm , 
Gpmax = 42,3 g, TaGpmax = 8000 h, Ag, = 0,004 kg/s [13] ) the values calculated by means of the 
relations (10) and (14) are the following: 
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Fig. 3 Values of the number k’ of the events N4 in function of considered engine load values 


163 


0,35 


0,3 4 
0,25 
=== 100% load 
-8-85% load 
0,2 4 
—t- 70% load 
€ 
a 
0,15 
0,1 
0,05 


0 FET T T — 1 
o © © © © © © © © © © © © © © © © © © © © © © 
O © © O © © © O O © © © O O O © O O © © © © 
AN O * nn © ~ on O m a M tN © p o O O m AN 
a — kam a m m" m" m" a a N N N 
t 
[h] 


Fig. 4 Occurrence probability of the number k’ of the events Nu, in function of considered engine load values 
and engine use duration period 


Summary 


As already stated, the worsening of operational characteristics of an engine which operates 
under partial load up to the instant t’ (Fig. 2), does not directly influence ship motion parameters. 

However such situation generates additional cost born by engine's user , which finally lowers 
his profit Z resulting from realization of the task. To determine potential losses due to increased 
fuel oil consumption is possible by assuming that for t = 0 the additional cost © = 0 and the profit 
Z takes its maximum valueZ(0)=Z,,. By determining the expected value and standard 


deviation of the random variable which describes the lowering of the profit Z (the increasing of 
financial losses) as: 


(15) 


the relation which describes the lowering of the profit Z along with time £ can be expressed as 
follows: 
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Z dla t=0 


x% 


Gpms — Sp Gyms — Sp 
= A A 
AU MECZY z r+o-Ag, at t| dla t>0 (16) 
AG AG 


Making use of the relation (16) one can determine , for a given instant t , costs generated as 
a result of increased fuel oil consumption; and the relation (14) makes it possible to determine 
occurrence probability of such number of the events F which will generate additional limitations 
during realization of the task (e.g. lack of possibility of arbitrary loading the engine within its 
working area), or will make its realization impossible at all. 
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Abstract 


In the paper has been make an attempt to replace the conventional system of fuel dose control and injection 
passing angle with elecrtronic control system, which has been realized this manner, that selected hydraulic 
accumulator, injector C-R type and fuel dose controller has been attached the the marine engine. 

It has been assumed, that by controlling the current impulse of controller it is possible to model whichever/any 
size of fuel dose, injected to the combustion chamber and to control the injection lasting time. 

The wave phenomenon and the flow loss in the fuel pipe of high pressure have been presented. 


Keywords: common rail, fuel pipe, wave phenomena. 


1. Introduction 


In the conventional and dividing injection pumps and in the fuel injection unit, the relation of 
pressing and fuel dosing process and the shaft or ring of bend rotation brings unfavorable change 
of the course of fuel injection parameter along with the change of the pumpshaft rotational speed. 

Broadening the scope of optimal engine operation becomes possible, thanks to applying the 
electronic control system, which enables controlling the course of basic dosing and fuel injection 
parameters in the whole scope of engine operation, in the various surrounding conditions, 
considering also fuel properties. 

The replacement of a conventional control system of fuel dose control and injection passing 
angle with the electronic control system in the test engine was realized this manner: 


** in the place of location the conventional injector was fixed selected in advance C-R injector 
produced by the Bosch Company 

** controlling the injection passing angle, size and multiplicity of the dose is realized thanks to the 
controller with the amplifier 

«œ fuel at demanded pressure is delivered by the pump through the hydraulic accumulator (Fig. 1). 


In the paper have been presented selected phenomena, which occur in the fuel pipe of high 
pressure at the moment of fuel injection. 


2. Scheme of the CR feeding system in the laboratory/test engine JSB type. 


The scheme of fuel system with attached common rail and fuel dose controller Fig 2.1. Similar 
feeding system was applied in the Sulzer RT — flex60C engine. By controlling the current impulse 
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of controller it is possible to obtain the whichever modelling of size and multiplicity of fuel dose, 
which flows to the combustion chamber in engine cylinder. 

The scheme of feeding system Common Rail in the one-cylinder test engine (JSB), which will 
replace the conventional fuel system, is presented below: 


Enaine operatina parameters 


Fig. 2.1. Scheme of the injection system CR of JSB engine: 1 — filter, 2 — low pressure pump, 3 — high pressure pump, 
4 — pressure sensor in the hydraulic accumulator, 5 — manometer, 6 — reversive valve with the opening control, 7 — 
overflow valve, 8 — signal amplifier, 9 — fuel dose controlling valve, 10 — pulverizer, 11 — low pressure fuel pipe 
(overflow), 12 — high pressure fuel pipe. 


3. Impact waves dispersion in the fuel pipes 


The wave flow in the centre, which fills the variable-section long pipe, is taken under 
consideration to describe the impact wave dispersion. The purpose is to explain the influence, 
which has the change of wave surface on the impact wave velocity. 

It was assumed, that the surface A(x) of pipe section slowly changes along its length (axis x) — 
little at a length in order of the pipe dimension, it is called hydraulic. It offers a possibility of 
applying appoximation (called hydraulic), so it can be assumed that all values in the stream are 
constant along every cross-section pipe, and velocity is directed along its axis, in other words, the 
flow is treated as a quasi-one-dimension. 

Such flow is defined by the equations: 


ov | ov 1p o 


3.1 
% Ge AG Fa) 
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OP ge (22) -0 (3.2) 
Ot Ox Ot Ox 
Op O 
A——+—(pvA)=0 3.3 
art (4) (3.3) 


The first of them is Euler equation, the second — adiabatic equation, and the third — continuity 
equation. 

To explain this problem, it is enough to consider the pipe, where the change of surface A(x) is 
not only slow, but absolute value of the change is low at the whole pipe length as well. Then, the 
stream disorder, related to the variable section, would be also low and equations (3.1) — (3.3) could 
be linearized. At last, the initial conditions should be imposed to eliminate appearing any foreign 
disorder, which could affect on the impact wave motion; only disorders related to the change A(x) 
are taken under consideration. The purpose will be accomplished/achieved, if it would be assumed, 
that impact wave initially moves at constant velocity through constant-section pipe and section 
surface changes only on at/on the right from certain point (which is assumed as x = 0). 

The linearized equations (3.1) — (3.3) have following form: 


0óv  Oóv 1 Oćp 
+V——+ = 
Ot Ox p Ox 
Oop | „0% z 0óp , „Odp -0 
Ot ox Ot Ox 
A 
Oóp p COR Só „BR -0 
Ot Ox ox A Ox 
where symbols without indexes mean constant values in the homogenous stream in the 
homogenous part of the pipe, and symbol 5 means a change of these values in the variable-section 


pipe. Multiplying the first and the third of these equations adequately by pa and a’ and summing 
up all three equations, can be written: 


a 8 wa” 064 
ae nee óv)=— nate 
(2 “(vale lb aa aay ae 


The general result of this equation is a sum of the general result of homogenous equation and 
the result of specific equation, which has non-zero right side. The first result is F(x — vt — at), 
where F is any function; it describes the sound disorder, approaching from the left side of the 
point, a — speed of sound. However, there are no disorders in the homogenous area, for x < 0, that 
is why it should be assumed that F = 0. So, the result amounts to the non-homogenous equation 
integral: 


0 


(3.4) 


2 
bm Tem 22: (3.5) 


The impact wave moves from a left to the right side at velocity v; > a; in the immobile centre 
with given values pı, pı. However, the motion of centre at the back of impact wave is defined by 
the result (3.5) in the whole pipe area at left side from the point, which discontinuity has reached at 
particular moment. After the wave crossing, all values in every of the pipe section remain constant 
in time, i.e. equal to the values, which they obtain at the moment of discontinuity crossing: 
pressure pa, density p2 and velocity vı — v2 (acoording to the symbols accepted in this paper, v2 
means velocity of gas in relation to the moving impact wave; its velocity in relation to the walls of 
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pipe is then equal to v; — v2). In these symbols (after dispering various terms of these values) 
identity could be formulated as follows: 


OA VI yt 
= Óp,+p Ovi— OV (3.6) 
{ Gis Eg ce zal 1 2) 


All values vı, v2, p2 could be defined by one of them, for IE óv,. To os end 
dependence variations could be described as 9,y;=,v2 and p,+p a= =p,+p 3v for the 


discontinuity (at given pı i pı) in the form of: 

P,Óvi = v20 P,+ PÔ v2» 2j(óv1—6w)= óp, + vð p, 
where j = p,y,j = p,vy» is an undisordered stream value. The following relation should be also 
included 


where derivate is calculated along the Hugoniot adiabatic curve. The calculations lead to a 
relation, which combine the change of impact wave velocity dv; in regard to the immobile gas 
before this change with the change of pipe section area, which means: 


1 OA _vi-v2tap 1+ 2v24, -h 3.7) 
A i Via? 1+h 
where the symbol is introduced again 
2 
n=- aP = pat: (3.8) 


2 
p 2 dp 2 dp 2 

The coefficient w, standing before the square bracket, is positive. Thus, the quotient symbol 

dv,/5A depends on the expression symbol in this brackets, for every stable impact wave is positive 


2dV2 
d 


2 


and when óv,/5A < 0. However, if any of discontinuity conditions j <-] and 


pee <j+2 2 , which are caused by corrugation, is fulfilled, the expression in brackets 
T dp, a2 
becomes negative and when 0v,;/5A > 0. 


Hypothetical shapes of the impact waves, 
moving along the pipe 


Fig. 3.1 
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The result offers possiblities of outlook explaining the origin of intability. Fig. 3.1 presents the 
“corrugated” surface of impact wave, moving to the right side; the directions of current lines are 
schematicly described by the pointers. During the impact wave crossing on the front protruding 
parts of the surface, the area 5A increases and on the back parts of the surface it decreases. In case, 
if óv;/6A < 0, it leads to the delay of the protruding parts and to the acceleration of the back parts 
and in relation to this the surface aims at smoothing. 

On the contrary, in case if óv,/5A > 0 the disorder of surface shape would demand: the 
protruding parts will be more protruding and the back parts will stay behind in higher extend. 


4. Fuel flow in the elements with constant axis-symmetrical-section area 


The d'Alembert solution, which was used by Allievi in the form of equations [5], was assumed 
to describe the unstable one-dimension fuel flow: 


plot) Po 1 Eey 6") f {1+ }er0(-n.€%) (4.1) 

a a a a 

u(x,t)=u,+ i r! A kóz)» i AGCA Kf) (4.2) 
ap a) ap a a 


a 


term, complemented byW. Bosch: 
x 
ep ke) (4.3) 


taking under consideration the flow loss with the momentum at length x. 


Symbols: 
£.{1-*)= Pi =] =apu,{t-*) (4.4) 
a a a 


- is an amplitude of pressure wave, which moves in the same direction as the flow at velocity a at 


length x; 
. 


- is an amplitude of pressure wave, which moves in the other direction than the flow at velocity a 
at length x; 


a=x/t (4.6) 
— speed of sound in the fuel; 
4T 
kx PA (4.7) 


— is a flow loss coefficient, considering the fuel mass momentum according to W. Bosch. 


P, — constant pressure before the disorder [N/m?]; 

u, - initial fuel velocity [m/s]; 

p -fuel density [kg/m*]; 

v - kinematic viscosity of fuel [m*/s|; 

A -flow section area[m?]; 

6 - fuel flow momentum coefficient (according to Burman and Deluca [1]). 
By introducing the same symbols as in the scheme, it is obtained: 

- in case ofl=x„-x=0 


P(xmt)=P„(t)=Po+ Pam(t)+ Pum (t) (4.8) 
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u(xmt)=um(t)=uotua()+ ttm (t) (4.9) 


-ifl=x,-x=0 
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Fig. 4.1. Scheme of unstable fuel flow in the axis-symmetrical pipe. 


P(x„t)=P,(t)=Po+ Pan(t)+ Pun (t) (4.10) 
U(X st) = un (t) = uo + Uan(t) + tm (t) (4.11) 
- where | = Xp - Xn > 0, 


I l 
PO=Prrem|- kuć | Pao #4 
a a 


+ ep km] . Pa + z) 


1 l 
EEEE: 
a a aj ap 


(4.12) 


a 


ap 
from the equations (8.60) and (8.62) follows relations: 
l l 
Palt)= erp kan}: Paf 1-7) (4.14) 
a a 
l l 
P„()=exp| kmó— |-Pm| t+- (4.15) 
a a 
l l 
Pin (t) = exp —kmm$ — "Pum ba (4.16) 
a a 


Volumetric stream of fuel, flowing through the sections m-m and n-n is defined by relations: 
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o, O= Apan O- Panl] (4.17) 
ap 
Q, (= = [Pan (t) Pan (e)] (4.18) 
ap 
V mn dP 
t)= t- a tt 4.19 
O,(¢)= Dy E- tm) =" (t) (4.19) 
Q,0--E|Pi(-tm)- Pan = tr) (0) (4.20) 
ap E dt 
where E is a fuel compressibility module [N/m?]. 
4.1. Fuel flow through the intermitten axis-symmertical-section of elements [2] 
ny rj 
Pall = L, Fit) 
CRON Q(t Qt) | at) alt 
uo ue ue A ub re 
P(t) | 55 pty à A, | w(t) 
5 A A i 
r mi > < > 
Fig. 4.2. Scheme of unstable fuel flow in the axis-symmetrical pipe. 
The courses of pressure waves in the appropriate sections are equal to: 
Pum (1) = Pam (1)+ Cr 1—Dnl|Pdm(t)}-Pum(t-ten b = Cn (4.21) 
Pan(t)= Pum (t)- C, {1-DmlP dm(t-tmn} Pun OF” +C (4.22) 
P.(0=P.()+C,U-DylPar(t)-Pus(-t,)" -C. (4.23) 
Pa (t) = Pis ()- Cs U-D,lPar(=trs)-Pus(O" T Cs (4.24) 
Substitutions: 
= Imi _ Ling 
bimn cm t ro (4.25) 
a a 
1 
=0.5¢° km 4.26 
Cn a P km ARE (4.26) 
k„+ 1 
= (03a p" — 4.27 
C, P (4.27) 
A (4.28) 
8 3k'-(K*Jk,)-2 | 
K+1 
C,=2a P (4.29) 
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Symbols: 

A,„A,„A,,A,„A; — surface area of appropriate sections [m?]; 

u,,U,,U, — mean velocity of fuel flow[m/s]. 

Fig. 4.3. The volumetric stream of fuel, flowing through the variable sections, is defined by 


equations: 
- throttles: 
ki (HA). J 
ch ee SARE | at, 
PO” BO 
Fig. 4.3. 
2 0.5 : 
0,0)-(2) (dusente (433) 

where: 


1 dla P,(t)> P,(t) 
senkl,0 dla p,(t)= P,(t) (4.34) 
-1 dla P;(t)< P,(t) 


(LA), -is an equivalent flow section between the sections k-1 [m?]; 
Qn Qn — are volumetric streams of fuel, which follows from the leak [m/s]. 


4.2. Mild change of direction — bend 


The Weisbach formula of loss coefficients for the circular-section bends is following: 


d 3,5 o 
€ =| 0,131+ (2) Q 
p 


90° 


In the table are given loss coefficient values, which are calculated according to above given 
formula for various angles a and relation d/p. 

The more detailed researches proves, that loss coefficient € depends not only on the relation 
d/p and the angle a, but also on the Reynolds number and the roughtness of the bend. 
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Fig. 4.4. Circular-section bend 


5. Conclusions 

Adapting the individual elements of electronic control system , which means an electronic 
control injector, a fuel dose controller with the amplifier and a common rail — for the conventional 
JSB engine there is a possibility of obtaining fully controlled process of electronic control. 

By controlling the current impulse it is possible to control size of fuel dose and also to model 
the manifold injection, which is characteristic of the electronic control injection systems. 

The simulation calculation were conducted basic of theoretical and empirical relations without 
verification and tests on actual model of valve and that is why they should be treated as a rated 
results. 

To verify the accepted model it would be necessary to perform empirical studies on the actual 
valve and use obtained results to revise the approved calculation model. 
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Abstract 


The paper deals with problem influence of changes variable stator vanes axial compressor settings of gas turbine 
engine on work parameters of compressor and engine. Incorrect operation of change setting system of variable vanes 
could make unstable work of compressor and engine. This paper presents theoretical analysis of situation described 
above and results of own research done on real engine. The next there are presented results of mathematical 
modelling of changes of gas turbine engine work parameters during change of angle setting of axial compressor 
variable stator vanes but in the most wide angle range than in real research. 


Keywords: gas turbine engine, axial compressor, variable stator vanes, modelling 
1. Introduction and purpose of research 


When in compressor construction is assembled system of setting change of variable stator 
vanes its task is made optimal cooperation engine units during permanent improvement of 
compressor characteristic. Perturbations in operation of this system could cause changes in work 
of compressor and engine similar like changes caused by changes of rotational speed or polluted 
interblades ducts of compressor. 

Compressor stage unitary work on radius is defined on base of equitation of angular 
momentum and has form (Fig. 1): 


lst = OF (Cu — Cry) = UAC, = UAW, (1) 


where: 
@— angular velocity, u — tangential velocity, r — rotor radius, 
Cium C2u — circumferential components of air stream absolute velocity on inlet and outlet rotor 
blades on radius r, 
Ac, Aw, — air stream whirl in rotor. 


That work is constant on whole depth of rotor blade. The sum of works is unitary work of 
stage. Involved change of variable stator vanes angle setting during kept at a constant level 
rotational velocity (constant u) caused change of air stream inlet angle in rotor vane /, (Fig. 1). 
It caused change of axial component of air stream absolute velocity on inlet cza what is equivalent 
with change of air mass flow m and change of air stream whirl Aw, in rotor. It influences 
on efficiency and work of stage [2]. 

Purpose of investigations made on real engine was determination influence of incorrect 
operation of axial compressor inlet guide variable stator vanes control system of gas turbine engine 
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on compressor and engine work parameters. 

Compressor characteristic is relationship between compression ratio 7's, compressor efficiency 
ns and air flow mass m and compressor rotational velocity n. It makes possible to determine the 
best condition of compressor and another engine units mating. The characteristic is using to select 
optimal conditions of air flow regulation and assessment of operational factors on compressor 
parameters. 

Therefore compressor should be so controlled in operational range of rotational velocity that 
the compressor and engine mating line has a stock of stable work. The main rule of compressor 
control during change of their rotational velocity or flow intensity is to keep up the stream inlet 
angles i value near zero. One of the most popular ways of axial compressor control is changing 
their flow duct geometry by application of inlet guide stator vanes or variable stator vanes of 
several first compressor stages [2]. 

This solution makes possible to change of air stream inlet angle on rotor blades of compressor 
stages by change of stator vanes setting angles during change of compressor rotational velocity. 
Fig. 1 illustrates, on example one stage of compression, rule of regulation of variable stator vanes. 


Fig. 1. Essence of control of compressor 's axial stage by changing the setting angle of stator vanes ring at 
changeable air flow velocity; a) decreased axial velocity, b) calculation axial velocity, 
c) increased axial velocity, d) schema of flow round of axial compressor rotor blades during constant rotor speed 
and constant air stream inlet angles; k — variable stator vanes ring, w — rotor vanes ring 


For average values with operational range of compressor rotor speed is situation on Fig. 1b — 
speed values and directions with subscript „Z“. In this situation is intermediate angle setting of 
stator vanes. Air stream inlet angle on rotor blades do not cause disturbance of stream flow by 
interblades ducts. For lower values of compressor rotor speed and in consequence lower values of 
absolute axial component velocity cy, , it is necessary to reduce the stream outlet angle of variable 
stator vanes a, (Fig. la). The angle reduction range should allow keeping the same value of stream 
inlet angle on rotor blades. Analogical situation takes place during work of compressor with higher 
rotational speed. For higher rotational speed absolute axial component speed cza”? increases. In this 
situation for keeping stable work of compressor and in consequence constant value of stream inlet 
angle on rotor blades, it is necessary to increase the value of stream outlet angle of variable stator 
vanes — Fig. lc. Application in gas turbine engine construction of control system of flow ducts 
geometry has a bearing on run inertance of unstable processes. 


2. Object of research 
The object of researches is type DR77 marine gas turbine engine, which is part of power 


transmition system of war ship. It is three-shaft engine with can-ring-type combustor chamber and 
reversible power turbine. 
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In compressor construction configuration of this engine there are used inlet guide stator vanes 
which make possibilities to change setting angle incidance (change of compressor flow duct 
geometry) in depend on engine load. This process is operated by control system which working 
medium is compressed air received from last stage of high pressure compressor. On Fig. 2 is 
presented block diagram of flow control signal of variable stator vanes system. 
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Fig. 2. Block diagram of stator vanes change setting mechanism LPC; CO — combustor, 
HPC — high pressure compressor, LPC — low pressure compressor, ŁK — variable stator vane, 
VIGV — variable inlet guide vanes 


Compressed air from last stage of high pressure compressor is supllied to working space of 
control actuator by cleaning and cooling block. Compressed air exerts pressure on control actuator 
elements. It causes moving of control piston which is connected with moving ring. This ring 
moves on circumference of low pressure compressor body. Ring is connected with variable stator 
vanes by levers. When the ring is moving stator vanes realize rotational motion changing the air 
stream outlet angle a, (Fig. 1). 

In cleaning and cooling block are holes. During researches air stream was bleeded by holes and 
less air was supplied to actuator. It caused change of setting angle a@xw of variable stator vanes. 
In consequence of that change flow duct geometry was changed. 

Experiment was carry out an engine load 0,5Pnom. For this load setting angle of variable vanes 
has value axw = - 4°. During change engine load in whole range from idle to full load setting angle 
axw of variable vanes changes in range from -18° to +18°. Realizing experiment a few parameters 
of engine work was measured and registered for three different setting angle akw of variable 
vanes: A— Qky = - 4°, B— agw = - 11°, C— agw = - 18°. Tab. 1 presents measured and registered 
in this same time parameters of engine work. 


Tab. 1. Parameters of engine DR77 work measured during researches 


Parameter Measurement range Unit Parameter name 
NPC 0 — 20000 [min ] low pressure rotor speed 
] 


high pressure rotor speed 

] power turbine rotor speed 

] subatmospheric pressure on compressor inlet 
] 

] 


Mie 0 — 22000 [min 
[ 


air pressure on low pressure compressor outlet 


[ 
[ 
[ air pressure on high pressure compressor outlet 


MPa] fuel pressure before injectors 


203 - 453 [K] air temperature on compressor inlet 
273 - 1273 [K] exhaust gases temperature on inlet power turbine 


3. Results of research 


Fig. 3 presents results of experiment. There are presented those parameters which are the most 
sensitive on change of LPC guide vanes setting angle. Change vanes setting from position A to 
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position C caused increase air flow resistance by stator vanes. In consequence of that 
subatmospheric pressure on compressor inlet p; decreases. It causes pressure decrease in next parts 
of compressor and engine flow duct (Fig. 3bc). In this way reduced air density flowing by 
compressor, for stable quantity of stream fule supllied to combustor, causes increase of 
compressors rotor speed. The most visible is increase of low pressure compressor rotor speed (Fig. 
3a) caused by directly influence on this compressor incorectly setting variable guide stator vanes. 


a) osit low pressure compressor rotor speed b) 03 - air pressure on low pressure outlet conpressor 
0,298 - 
'g400 z 
= 
5 2 0,296 + 
9300 4 a 
o a 
a 0,294 - 
s 

9200 - 

0,292 - 

9100 TT 0,29 $44 +, 
-20 -18 -16 -14 -12 -10 -8 -6 -4 -2 0 -20 -18 -16 -14 -12 -10 8 6 4 2 0 
angle setting of controllable vanes æ xy [°] angle setting of controllable vanes agy [°] 

c) air pressure on high pressure outlet conpressor 
0,85 T Tol a a 
I nl zb KAY al 
b+ sel tate Paci zj” z 
I nl [ofl atte] 
-+ +4 +1--+1 
I nl © deg dd 
2 (hie AB 3 
- + +4 +1--+4 
Sos ENGL 
a 
a 
0,75 + 


-20 -18 -16 -14 -12 -10 -8 -6 -4 -2 0 


angle setting ofcontrollable vanes gy, [*] 


Fig. 3. Change of engine work parameters values in function of variable inlet guide stator vanes setting angle: 
A- agw =-4,B-Qgy =-1I', C- agw =- 18° 


Gasodynamical connection between low pressure compressor and high pressure compressor 
absorbs disturbances work of low pressure compressor which are transferred on high pressure 
compressor. Therefore range of change high pressure compressor rotor speed is lower than low 
pressure compressor. In this experimental it is below 1% and it is in measuring error range. 

Change of subatmospheric pressure is above 5% stable value of this parameter. Changes of low 
and high pressure compressor outlet presure are adecuately above 1,3% and above 2,4% 
undisturbed value of angle setting akw = - 4°. 

Changes of pressure and air mass flow intensity values accompanied disturbing work of 
compressor, during constant fuel mas flow intensity in combustor, caused enrichment of fuel 
mixture. As a result of that, temperature combustor outlet gases increases. In experiment was 
confirmed tendency changes of gases tempertaure values even though range of thoses changes is in 
measuring error range. 

On the base of results of experiment there were determined mathematical equations modelling 
changes of particular engine work parameters in function of variable inlet guide stator vanes 
setting angle ay: 


180 


(2) 
(3) 
(4) 
(5) 
(6) 
(7) 


to +18” were 


0,16330axkw + 526,33 


0,74490kw + 2,6020w + 9234,5 
0,0204axw” - 1,1224axw + 12598 
pi =-10°axw? - 10° ay + 0,0077 
10 any? + 0,0029 xy + 2,9814 


pa = 2:10 aw? + 0,01430xw + 8,1771 


0,0204axKw? + 


In this case range of change of variable inlet guide stator vanes 
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Fig. 4 presents results of solution of an mathematical equations describing of changes of engine 
work parameters values. Modelling was cary out an state engine load what was equivalent 
not possilble to realize on real engine. It is caused by technical restrictions on the engine. 


setting angle axw was widen from -18° to +18°. Researches in range aw from -4° 


unchangable fuel mass flow. 


gotten during mathematical simulation 
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Fig. 4. Change of engine work parameters values in function of variable inlet guide stator vanes setting angle 


4. Conclusions 


Change of variable inlet guide stator vanes setting angle a@xw from -4° to +18” caused increase 
of flowing air stream outlet angle a, (Fig. 1). It decreases air flow drag on low pressure 
compressor inlet that caused decrease of subatmospheric pressure. During keeping constant engine 
load (constant fuel mass flow) absolute axial component velocity c; increases. It exerts an 
influence on air mass flow m increase. Simultaneously the absolute axial component velocity cja 
increase caused decrease of air stream whirl in rotor Aw,. The effect of above is reduction of the 
compressor stage unitary work — equation (1). In consequence of that low pressure compressor 
rotor speed increases (Fig. 4a). In connection with decrease of subatmospheric pressure it caused 
increase of air pressure on low pressure outlet compressor (Fig. 4d). In spite of slight decrease of 
high pressure compressor rotor speed the increase of air pressure on low pressure outlet 
compressor involves increases of air pressure on high pressure outlet compressor (Fig. 4e). This 
slight decrease of high pressure compressor rotor speed caused increase of gases flow drag in next 
gas turbine engine units for the combustor. The effect of above is slight increase of exhaust gas 
temperature on power turbine inlet. 

Multi-shaft construction of gas turbine engine reduces effects of incorrectly setting of variable 
stator vanes. Therefore compressors of three-shaft gas turbine engine do not require variable 
stators vanes as many stages as compressor of two-shaft engine with the same achievements. 

Preliminary researches confirm necessity for making inspection of correct operation of variable 
stator vanes system control. It makes possibility of elimination this factor from group of factors 
informing about deteriorating technical state of engine which are identified during diagnostic 
inspections. 
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Abstract 


The article provides the specification of the experimental testing method of the heat exchange between the 
fluidized bed and the heating surface immersed therein during ship s movements on the waves causing interference. 
The research works have been conducted on physical model of fluidized bed boiler positioned in the cradle simulating 
the vessel’s motion on regular wave of sinusoidal type. The design features of the test stand make it possible to 
conduct the testing for the processes occurring in the boilers both with bubbling and the circulating fluidized beds. 
The method of determination of the heat transfer coefficient has been specified. The presented research results of the 
heat transfer coefficient have proven that the conditions of heat transfer are of explicitly local nature strongly related 
with the mass of the material present in the fluidizing column, flow-through conditions, fluidizing column inclination 
angle from vertical and swinging movement of the platform. The testing results encourage to continue the research 
works on a larger scale aiming to verify the hypothesis that in the shipboard conditions the circulating fluidized bed is 
less susceptible to the ship’s rolling than the bubbling fluidized bed and that the significant conditions of heat 
exchange are not impaired. 


Keywords: environment protection, heat exchange, ship, fluidized bed boiler, coal 


1. Introduction 


Diesel engines prevail in the modern shipboard propulsion systems. This predominant position 
results from their high efficiency which in case of two-stroke low-speed engines exceeds 50% as 
well as the possibility of burning low quality cheaper fuels. However, while developing the future 
ship’s propulsion systems one needs to account for the decreasing oil resources and the experience 
proving the strong influence of some political events or natural disasters on the oil prices. Thus the 
research in the field of application of the alternative fuels is of urgent nature at present. The basic 
role here is played by coal which already in the past fuel crises proved to be a sound alternative 
fuel. Its resources are sumptuous as according to the USA Energy Information Administration with 
the current global energy consumption maintained the coal resources would suffice for 600 years. 
Also in terms of price the present cost of energy production from coal is approximately 6 times 
less than that from oil [6]. 

However, a serious deficiency of coal as a fuel is high emission of CO2. As much as 41% of 
global emission of this compound in the energy production comes from coal combustion. Within 
the land power engineering it is planned therefore to apply systems separating CO» to store it in the 
earth and to apply advanced power engineering technologies. The fluidizing reactors within IGCC 
are considered as particularly attractive in view of their high efficiency and purity of the energy 
production [5]. Because of the size of the installations and degree of complexity of the system, it is 
not suitable at the present stage to be applied on board ships. However, coal application can be 
assumed as fuel on ships by its combustion in high efficiency fluid-bed boilers. The problem of 
CO» emission may be solved on the other hand in the future by its separation from exhaust gases 
and subsequent release into the water [4]. The studies of the possibility of application of these 
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boilers on ships have been assumed by the end of the 1980’ of the last century in Germany [3]. In 
view of the shortage of the information on the experimental research in this respect a test stand for 
the pilot testing of the physical model of the “shipboard” fluidized bed boiler has been built in the 
Department of Heat Engines and Marine Power Plants of the Faculty of Maritime Technology of 
the Szczecin University of Technology [1]. 


2. Measurement Apparatus and Technique 


The specially designed and executed test stand enables the testing to be conducted for the 
processes occurring in the boilers both with the bubbling and circulating fluidized beds in the 
conditions simulating the ship’s motion on waves. The detailed diagram and the specification of 
the stand has been provided by the author in the study [1]. The Fig. 1 shows the view of the stand 
while the column is inclined (a) and the interior of the fluidizing chamber (b). 


=~ 
3 fluidizing 
A column 


thermocouples 


heating 
element 


bed 


cradle 


actuator 


a) b) 
Fig.1. View of the stand 


A static method has been adopted for the investigation of the local heat transfer coefficients 
between the fluidized bed and the heating surface immersed therein. This method is based on the 
energy balance in the stabilised condition of the set-up: heating surface — fluidized bed. A pipe 
with electric spiral (electric heater) forms the heating surface. The determined value of the local 
coefficient of heat transfer aay is a mean value that can be obtained from the Newton's equation 

Q 
Oy ==, (1) 
where: 

Ò — heat stream supplied to the set-up, 

A — heating pipe surface, 

At — difference of the mean temperatures of heating pipe surface and fluidized bed. 


The heat stream is determined on the basis of the measurement of the electric current supplying 
the heating spiral placed inside the pipe. In the preliminary investigations the provision has been 
made for the temperature readout in 6 points of the bed above the heater, the points located in two 
rows in the plane of symmetry of column and heater. In the row situated closer to the heater three 
thermocouples are positioned, and in the farther row three remaining ones. The temperature of the 


184 


heater wall is measured by use of thermocouple made fast to its surface from the top. As the pipe 
is made of copper thus the material of high heat conductivity coefficient, a uniform temperature 
distribution all over its surface has been assumed. Additionally the interior of the pipe where the 
electric spiral is placed has been filled with sand thus increasing its thermal inertia and the uniform 
heating of the wall. Teflon bar/plugs have been inserted in both ends of the pipe to reduce the heat 
losses from the fore walls. Heat transfer coefficients have been calculated as average values from 
at least three temperature readouts upon reaching the stable condition by the set-up. The layout of 
heater and thermocouples inside the column which are also visible in the Fig. lb is shown in 
Fig. 2. 
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Fig. 2. The positioning of heater and thermocouples in fluidizing column 


Poppy seeds have been used for the purpose of the pilot testing of heat transfer coefficient at 
the stand in question. The seeds have the shape similar to ball of diameter d = 1 mm and the 
density nearly three times less than sand, ie ps = 714 g/m’. The velocity of the beginning of the 
fluidizing (Vikr) for poppy seeds is 0.17 m/s. Owing to this it has been possible to apply lower 
capacity blower and eliminate the risk of scratching the column. 


3. The Results of Testing and their Analysis 


The described method has been used at the erected stand to conduct several pilot measurements 
in respect to various aspects of bed behaviour. These allowed to evaluate the model suitability for 
conducting the experiments on a larger scale in the future. The first test aimed at determination of 
the local heat transfer coefficients for the bed of stationary height H=12 cm , in the function of 
fluidizing air jet velocity. The bed transition into fluidized state has been observed at air velocity 
of v;=2.2 m/s. The air velocity has been increased in steps until the moment when the grains were 
raising from the column which was observed at the velocity of 6.2 m/s. Since that moment the 
transition to the circulating fluidized bed took place owing to the system of grain reverse in the 
model. The heat transfer coefficients have been determined on the basis of the measurements for 
each velocity, upon reaching the stabilised condition by the set-up. The power of the heater during 
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the experiment has been maintained at the constant level of O= 45 W. Figure 3 presents the values 
of the local heat coefficients in the column axis at the height of 25 and 55 mm above the heater 
(points 2 and 4 in diagram in Fig. 2). 

Initially the coefficient value increase followed by the marked decrease has been observed 
along with the increase of the air velocity. The moment of the rapid drop in the values of the heat 
transfer coefficients occurred as the deposited material arose and began to circulate. This is 
explained by the fact that air velocity increase is accompanied by a large decrease of material 
concentration [2]. As evident from the graph of the local values of the heat transfer coefficients in 
points 2 and 4 differ slightly. Still somewhat bigger value of the coefficient in the point closer to 
the heater prevails 

For the comparison the same figure shows the values of the local heat transfer coefficients in 
clean air without poppy seeds (H=0). As can be seen the heat transfer coefficients in the bubbling 
phase are significantly higher and only upon its dilution and transition into the circulation phase 
these get closer to the values achieved in clean air. The heat transfer coefficients in air grow in the 
obvious manner together with velocity, the value of the coefficient being somewhat higher at the 
point closer to the heater. 
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Fig.3. Local heat transfer coefficients for the fluidized bed and air in the function of air velocity 


A similar experiment has been conducted for H=9 cm in order to establish the influence of the 
ratio of the stationary bed height to the column diameter (H/D) so as to allow in the future to 
record the test results pursuant to the theory of similarity by means of the dimensionless equation, 
eg of the following type: 


Nu=fRe~), 2) 


where: 
Nu — Nusselt’s number, 
Re — Reynold’s number. 


The results of the measurements of the local heat coefficients for the beds of stationary heights 
H=9 cm and H=12 cm are shown in Fig. 4. 
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Fig. 4. Local heat transfer coefficients for fluidized bed in the function of the air velocity for poppy seeds for the 
stationary bed heights of H=9 cm and H=12 cm. 


It shows that for bigger stationary height of the bed in the same points higher values of local 
heat transfer coefficients occur. This can be explained by the higher poppy seed concentration in 
way of the measurement points. Because of the scarce number of experiments done so far with 
various stationary heights this tendency cannot be regarded as a lasting or permanent one. With 
bed heights even bigger the fluidizing nature is likely to be altered, e.g. channelling would occur 
which would influence the conditions of heat exchange. 

The investigation of the fluidized bed during the constant inclination of the column has been an 
experiment related to the conditions likely to occur in the boiler operation on ship. Figure 5 
presents the hydrodynamic characteristics of 9 cm high (in vertical) stationary bed for the constant 
inclination angle of 30°. In these conditions the bed gets fluidized faster. This is caused by the out- 
of-parallel location of the bed surface in relation to the grid. The bed movements have been 
observed already at the velocity of approximately 1.4 m/s in its “more shallow” part. A distinctive 
circulation of the grains has been noticed as the fluidizing air velocity has been increased upon the 
bed reaching the full bubbling fluidized condition. The produced circulating movement is of the 
dextrorotary type (right-sided) with the column inclined to the right and just the opposite — 
levulorotary type (left-sided) with the inclination to the left. Figure 6 demonstrates the movement 
of the material in column during inclination to the right. 
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Fig. 5. Hydrodynamic characteristics of the bed (H=9cm) with the column inclined by 30° 
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1 — separating grid 
1 —7 2 — bed surface in stationary condition 
R 3 — bed material circulation 
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Fig. 6. Diagram of the grain movement in bed 


For the purpose of comparison the investigation of the local heat transfer coefficients has been 
conducted for the column both position in vertical as well as inclined for the bubbling bed. The 
results of these tests are shown in Fig. 7. 
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Fig. 7. Local heat transfer coefficients for the fluidized bed in the function of velocity for poppy seeds with the column 
positioned in vertical and inclined by 30° 


The experiment indicates a marked decrease of the local heat transfer coefficients in all 
measurement points for the column inclined. For the sake of clarity the coefficient values are 
presented in the axis of column, i.e. in points 2 and 4. The decrease in the value of the coefficients 
is caused most probably by the poppy seeds circulating movement which has been produced in the 
bed. 

The figures 8, 9 and 10 present on the other hand the values of the local heat transfer 
coefficients determined along the parallel line to heater axis at the distance of 25 mm above it 
(measuring points 1, 2 and 3) for the vertical column (Fig. 8), column inclined by 30° to the right 
(Fig. 9) and inclined by 30° to the left (Fig. 10). In every case values are presented for two 
different air velocities. The mass of the loose material in both cases has been the same and 
corresponded to the height of stationary bed in the vertical column H=12 cm. In all cases it can be 
seen that the biggest value is reached by the heat transfer coefficients in the column axis thus over 
the heater centre. Slightly less values achieved at the walls may result from the circulation inside 
bed caused by wall-related effect. 

Comparison of the values of the heat transfer coefficients for three column positions looks 
interesting. Column inclination from the vertical results in reduction of the value of the heat 
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transfer coefficient at the side towards which the inclination took place. This may be the effect of 
the differentiation of poppy seeds concentration on both sides of the inclined column. 
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Fig. 8. Values of the local heat transfer coefficients at the line of 25 mm above heater for the vertically positioned 
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Fig. 9. The values of the local heat transfer coefficients along the line 25 mm above the heater for the column inclined 
by 30° to the right 
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Fig. 10. The values of the local heat transfer coefficients along the line 25 mm above the heater for the column 
inclined by 30° to the left 
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The last testing of the bed the results of which are illustrated in Fig. 11 consisted in the 
determination of the mean values of the local heat transfer coefficients at the parallel line to heater 
axis 25 mm above it (measuring points 1, 2 and 3) during constant pendular motion of the column 
in the period of 28 seconds and with the air velocity 3.5 m/s. Getting the column into pendular 
motion caused the decrease in the mean values of heat transfer coefficients in these points. This is 
explained by the temporary cyclical dilution of the bed in accordance with the column swing 
period. 
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Fig. 11. The values of the local heat transfer coefficients along the line 25 mm above the heater for the column 
constantly positioned vertically (alpha V) and during constant pendular movement with the 28 s period (alpha P) 


4. Conclusions 


The presented results of the pilot testing confirm the suitability of the model and the applied 
method for the bigger scale research. The conditions of the hear transfer are of visibly local nature 
depending on the mass of the loose material present in the column, on flow-through conditions, on 
the angle of column positioning and its pendular movement. Further investigations are necessary 
to determine the relation for the specification of the heat transfer coefficient. From the results 
obtained for the vertically positioned column which are confirmed in literature it is deemed 
purposeful to conduct research works mainly for various column positions as well as for the 
column being in constant motion both with bubbling and circulating fluidized bed. 
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Abstract 


The results of the investigations of exploitation emissions the harmful components of the fumes of diesel engines were 
introduced in the work. Obtained results were subjected to a statistical study according to new computer procedures. 
Qualitative and quantitative reports were established for the level and kind of emission in reference to the changes of 
the state of studied engines. 


Keywords: combustion engines, toxic fume components, exploitation investigations, environmental protection 
1. Introduction 


The assumption for researches of this work was the performance of the analysis influence of 
engine warming on the harmful emission at these states of engine’s work, especially concerning 
climate conditions in Poland. In the range of researches, the analysis of the emission of harmful 
compounds was performed during the warmed up engine in the neutral race at different 
temperatures of the environment. 

The conducted considerations indicate that researches of the work of high-pressure engines of 
different destination in the conditions of exploitation take place usually in changeable unsettled 
conditions, or are a majority of its working time, which considerably influences the general 
emission of harmful components of fumes. 

The contribution of harmful components of fumes diesel engines into total atmosphere 
pollution is as follows: there are mostly solid particles (PM) and nitro oxides (NOx) in fumes, 
whilst in smaller amounts there is carbon oxide (CO) and not burned hydrocarbons (HC). 

The results of realized in exploitation tests on a chosen group of diesel engines allow to 
determine, practically and cognitively, important premises in the field of toxic effects of diesel 
engines on the environment. 


2. Research objects 


The researches of his work, in the field of recognizing toxic components generated by 
diesel engines for different technical states and change able external temperature, were performed 
on a engine S-359. Results of these researches in the range of quality changes of the toxic 
components of fumes in operation conditions was performer during 2 years of exploitation on the 
group of 20 vehicles with such engines. The units of the investigative position were showed on 
Fig. 1. 
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Fig. I General view of test stations 


The object of research in this work was S-359 engine with self-acting fusion whose basic 
technical data is presented in Table 1. It is an engine of a wide practical application, and 
characterized by small unitary fuel use, good dynamic characteristics and small damageability. 


Tab.1. Basic technical data of the engine S- 359 [105] 


Cylinder formation row, vertical 
Number of cylinders 6 
Cylinder diameter 110 mm 
Piston stroke 120 mm 
Swept capacity 6,842 dm” 
Compression degree 17 
Order of cylinder work 1-5-3-6-2-4 
Maximum Power 110 KW with 2800 min” 
Maximum turning moment 438Nm with 1800-2100min' 
Minimum unitary fuel use 224 g/kWh 
Statistical angle of pumping beginning 18,5°OWK before GMP 
Injection system Direct 
Injection pump P-76G10 
Injection pressure 22MPa 


The engine is a running unit for trucks: Star 200 — street, Star 266 — cross-country, 
produced in Factory in Starachowice (at present: Star Trucks Sp.z 0.0). These cars are widely used 
in the national industry, as well as military service. 

The tested combustion engines belong to the group of exploitation objects, used in difficult 
training conditions of military service. Large and changeable loads of engines implied by 
inexperienced drivers diversified their technical state, which for the researches of his work posed a 
challenge in the range of preparing the experiment, its proper realization, and careful concluding 
and statistical work. 


3. The conditions of exploational investigations 


It mattered considering the acquisition of different temperatures, in which the engine S-359 
was thermally stabilized, and considering the temperature of the air used for running the engine. 
Before proceeding with the tests, the following were checked and regulated: 

a. technical state of the engine, 

b. injection pump at the probing station type PW-8, predestined for testing fuel equipment of 
high-pressure engines with regard to dosage, performing, according to BN-88/1301-16 
velocity characteristics of fuel injection, 

c. injectors used for the tests were checked and regulated on an injector probe type PRW-3, 
performing the evaluation of pressure of the injector’s opening, tightness and trickling of 
the sprayer, and the correctness of fuel spraying, 
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suction and exhaustion valves — according to the manufacturer’s suggestions. 

During the test, the following were registered: 

multicomponent composition of exhausted fumes of the engine, 

smoking of fumes with a smokemeter AVL. 

Fume tests with respect to the quantity of toxic substances were performer with the use of a 
multicomponent analyzer of fumes LANCOM, whose general image is presented on the Fig. 2. 
The analyzer LANCOM enables the measurements of: CO, CH, NOx, SO», fumes temperature and 
environment temperature. 
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Fig. 2. General image of fume analyzer LANCOM with fume acquisition probe 


The measurements of smoking degree of fumes diesel engines were performer with the use of a 
smokemeter AVL-4000 for the need of further statistical processing, the measured values were 
saved in a sheet (Excel). Exploitation tests were performed in real conditions on the group of 20 
vehicles. 


3.1 Testing conditions 


The engine, before each test, was subjected to thermal stabilization, thanks to which all 
elements of the engine and exploitation liquids and exhaust system had the same temperature (for 
every test) equal to the environment temperature. Environment temperature and motor oil 
temperature were measured directly prior to each measurement, and if the temperature differences 
did not exceed 1°C, the measurement began. Also performed were tests in the conditions of a hot 
start-up. i.e. during a start-up of an engine beforehand warmed up to a normal temperature (oil 
temp. 80°C) in certain environment conditions. During the measurement, registered were 
(LANCOM, AVL) the contents of carbon oxides (CO), hydrocarbons (HC), nitro oxides (NOx) in 
fumes, motor oil temperature, rotational speed of the crankshaft, environment temperature, and 
fumes smoking. 

Exploitation tests were realized in real conditions on a group of 20 randomly chosen 
vehicles with the engine S-359 throughout the period of 2 years in a military facility. The tests on 
toxic substances emission in fumes (CO, HC, NO, and PM) of diesel engines during a hot start-up 
(oil and cooling liquid temperature — 70-80°C) were performer at different environment 
temperatures, with a special consideration of run kilometers of vehicles, and operation- 
maintenance actions carried out during breaks in tests. 

Average daily run of each car was about 70km, which with sustaining the regime of car 
operation suggested by the manufacturer gives on average about 14000-16000 km per car yearly. 
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3.2 Exploitation tests results 


From the general number of 20 cars subjected to exploitation tests, 10 cars were chosen for the 
initial analysis, for which checked was the effect of the quantity of emissions of toxic fume 
components as well as the degree of smoking in relation to run kilometers in the testing period. 

The chosen vehicles (num.4) were evaluated in respect of the contents of CO, HC, NOx in 
fumes, and smoking of the hot engine, which is shown as an example in the Fig. 3. 
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Fig. 3. Overall values of CO, HC, NO, emissions engine 4 in relation to run km 


The following Fig. 4 presents the juxtaposition of values of engine (vehicle num.4) fumes 
smoking for different km runs of the car. 
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Fig.4. Smoking of 4 engine fumes depending on the km run 


The vehicle (num.10) were evaluated in respect of the contents of CO, HC, NO, in fumes, and 
smoking of the hot engine, which is shown on the Fig. 5. 
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Fig. 5. Overall values of CO, HC, NO, emissions engine 10 in relation to run km 
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The following Fig. 6 presents the juxtaposition of values of engine (vehicle num. 10) fumes 
smoking for different km runs of the car. 
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Fig. 6. Smoking of 10 engine fumes depending on the km run 


Even brie analysis of the presented results of these tests indicates a visible increase of toxic 
fume components and engine smoking together with increasing mileage of the car. Other test 
results of the whole group of 20 cars confirm it visibly. It can be stated that together with the 
increase of the number of run kilometers, the quantities of carbon oxide (CO), hydrocarbons (HC), 
nitro oxides (NOx), and smoking go up. 


3.3 Results investigations regress for the exploational investigations 


Analysis of the regress for results got in exploational investigations 10 cars were conducted 
stay for selected about the largest course. Put given these in the Tab. 2 and Fig. 7 represents the 
results of polynominal regressions, together with coefficient R”. 


Tab. 2 The results of exploitation investigations for chosen 10 vehicles 


ARAMETR 
Co HC Nox k 
NR VEHICLES 
1 0,567 0,797 0,978 1,010 
4 0,141 0,135 1,000 0,945 
5 0,096 0,151 0,389 1,000 
8 0,053 0,101 0,153 1,000 
9 0,064 0,132 0,101 1,000 
10 0,311 0,638 1,000 0,989 
11 0,169 0,191 0,311 1,000 
15 0,021 0,056 0,044 1,000 
18 0,059 0,156 0,250 1,000 
19 0,173 0,195 0,461 1,000 


y = -0,0065x + 0,0567x- 0,0292 y = -0,0083x" + 0,0719x-0,0011 y= 0,004x - 0,1647x + 1,0558 
R? = 0,8007 R? = 0,1375 R? = 0,7631 
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Fig. 7. The date for the polynominal regressions analysis of the studied group of cars 
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4. Fume test results for engine run on BIO-D10 


The possibilities to satisfy increasingly strict regulations force vehicle manufacturers to 
search new solutions. One of the ways is looking for new, ecologically purer, fuels, among which 
galenic fuels begin to play a dominant role. 

Vegetable oils in their pure form, also colza oil, are not suitable for engines with self- 
exciting fuse, mainly because of their increased density and viscosity, low cetane number, and 
insufficient immunity to low temperatures. These disadvantages are absent in products of chemical 
processes of vegetable oils, called methyl esters, which, combined with diesel oil in appropriate 
proportions, are called Bio-diesel. 

Taking the above into account, quality comparison between diesel oil (ON) and the oil 
BIO-D10 was carried out during test realization. The comparison of qualities of BIO-D10 fuels 
with ON, used for tests, is shown in the Table 3. The characteristics of BIO-D10 fuel, produced in 
refinery Trzebinia are included in the certificate No. 5100634 from 7.10.2005, issued by the 
manufacturer. 


Tab. 3. Comparison of chosen characteristics of test fuels of the engine 5-359 


Characteristics Unit Diesel Oil BIO-D10 
Density at temp. 15°C g/m 0,836 0,841 
Kinematical viscosity at mm/s 2,76 2,82 
temp. 40°C 
Fuse temperature °C 63,5 72 
Cetane number 51,1 52,4 
Sulfur contents mg/kg 6,9 4,8 
Temperature of cold filter -30 -25 
blockage °C 
Remains after %(m/m) 0,002 0,003 
incineration 
Water contents mg/kg 68 93 


Juxtaposition for easier comparison of emissions of toxic fume components of diesel oil 
and BIO-D10 during cold and hot start-ups of the engine in diversified environment temperature is 
presented in the following figures. 


Juxtaposition of CO emissions values of the engine S-359 
run on ON and BIO-D10 
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Fig. 8. Juxtaposition of CO emissions values of the engine S-359 run on ON and BIO-D10 
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Juxtaposition of HC emissions contents of the engine S-359 
run on ON and BIO-D10 
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Fig. 9. Juxtaposition of HC emissions contents of the engine S-359 run on ON and BIO-D10 


Juxtaposition of NOx emissions contents of the engine S-359 run 
on ON and BIO-D10 
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Fig. 10. Juxtaposition of NO, emissions contents of the engine S-359 run on ON and BIO-D10 


Juxtaposition of smoking of the engine S-359 run on ON 
and BIO-D10 
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Fig. 11. Juxtaposition of smoking of the engine S-359 run on ON and BIO-D10 


As the result of performed tests and result analysis, it was concluded that using the fuel 
BIO-D10 instead of diesel oil at cold engine start-up, at environment temperature of 5°C, causes 
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the decrease of emissions of CO with 16%, HC with 13% and smoking with 6%. Whilst the 
emission of NOx grows with 10%. 

It can also be stated that with the use of BIO-D10 at cold and hot start-ups at environment 
temperatures of 10°C and 20°C, the contents of CO, HC and smoking in fumes decreases. Whilst 
the contents of NOx go up. 

The presence of oxygen in the fuel causes the increase of nitro oxides with simultaneous 
decrease of carbon oxide and hydrocarbons, thus giving easier possibility to regulate toxicity of 
fumes by delaying the angle of fuel injection initiation. Steering the combustion process can be 
performer in a much wider range. 


5. Summary 


The need for detailed analysis of phenomena changes of state destruction of the examined 
engines, with a vast number of measuring data, requires the use of specialized methods of 
statistical concluding. The presented results were submitted to statistical analysis, where the 
methods OPTIMUM and AVD were used, as well as correlation and regression methods. It gave 
the possibility of quality and quantity comparison of results of fumes contents from stationary tests 
and exploitation researches. The results of this research allow a model (mathematical relations) 
determination of relations between smoking and the quantity of toxic fume components of a high- 
pressure engine. 

The performed tests and analyses in his work’s researches indicate to the conclusions: 

1. In the engine of self-acting fuse (ZS), the emission of carbon oxide (CO), hydrocarbons (HC) 
and smoking are considerable, especially during start-up and engine warming. 

2. Along with the decrease of environment temperature, the emission of CO, HC and smoking 
increase, whilst the quantity of NOx goes down, providing premises confirming the specified 
regulations of forming dangers on the side of engine fumes emission. 

3. The phases of start-up and warming up of the ZS engine are characterized by increased fuel 
usage and increased emission of carbon oxide — CO, giving information and sensitizing vehicle 
users to these harmful for the engine working conditions. 

4. The influence of environment temperature on the emission and smoking of fumes during hot 
start-ups is weaker than during cold start-ups. 
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